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Abstract PHOTOS Monte Carlo is widely used for simulating QED effects in decay of intermediate particles 
and resonances. It can be easily connected to other main process generators. In this paper we consider decaying 
processes 7* — > TT + TY~( r y) and K ± — > ir + ir~e ± v (7) in the framework of Scalar QED. These two processes are 
interesting not only for the technical aspect of PHOTOS Monte Carlo, but also for precision measurement of 
o,qbd{Mz), 3 — 2, as well as rnr scattering lengths. 
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1 Introduction 

In high energy experiments, one of the crucial 
works is to compare new experiment results with pre- 
dictions from the theory. If the agreement is obtained, 
the theory is proved to be true. Otherwise one may 
think that the theory calculations turned out to be 
wrong or the effect of new physics appeared. Monte 
Carlo generators, rather than analytical calculations, 
are required to provide theoretical results of real ex- 
periment interest. The PHOTOS Monte Carlo 0, Q is 
a universal Monte Carlo algorithm that is designed 
for simulating QED radiative corrections in cascade 
decays. It is widely used in low energy experiments 
and high energy experiments as well. The program 
is based on exact multiphoton phase space while the 
matrix element is approximately taken as process in- 
dependent multidimensional kernel. In some cases, 
the exact first order matrix element is employed in 
order to improve the precision. 

From the experience gained in KKMC project 
0, 0] where spin amplitudes were used, one concludes 
that spin amplitudes are essential for design and tests 
of the program, in particular for choice of the single 
emission kernels. The analysis of the spin amplitudes 
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and tests for the algorithm in case of Z decay into 
pair of charged fermion was given in Ref. fl. The 
case of the scalar particle decay into pair of fermions 
is covered in Ref. @ and the decay of spinless parti- 
cle into pair of scalars is studied in Ref. [7J . The case 
of W decay was covered in Ref. @, Q . 

In this paper we will study spin amplitudes for 
7* — > 7r + 7r~(7) decay. It not only provides example 
for studies of Lorentz and gauge group properties of 
spin amplitudes and cross sections, but also is impor- 
tant to improve theoretical uncertainty of PHOTOS for 
this decay. Furthermore, it is of great experiment in- 
terest since its relevance to precision measurement of 
oc Q ed{M z ) and g -2. 

K e4 decay could give the unique information on 
the value of s— and p— wave tttt scattering lengths. 
The high statistics measurements of K e4 decay has 
been performed by NA48/2 collaboration at CERN 
[loj . Theoretical predictions are well calculated 11 1 
and the difference between theoretical prediction and 
experiment value is about 7%. In particular, QED 
corrections to this process are known to be non- 
negligible. They need to be taken into account with 
the help of Monte Carlo because their size depend 
on detector acceptance. In NA48 experiment, to take 
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into account QED effects, PHOTOS Monte Cario is used 
together with Couiomb correction (see Ref. 10(). In 
this paper we consider QED radiative corrections to 
K e4 decay and implement it into PHOTOS Monte Carlo. 

The paper is organized as follows: process 7* — > 
7r + 7r - (7) is studied in Section 2, Section 3 will come 
to K e4 decay, summary will be in Section 4. 

2 7* — > n + 7v~ (7) 

In order to match the parton shower and the 
hard bremsstrahlung matrix element into interference 
weight, the matrix element of 7* — > 7r + 7r~ (7) needs to 
be studied in great detail. Its gauge invariant parts 
need to be identified and relations to amplitudes of 
lower orders have to be found. 

2.1 Spin amplitude 

Consider the process e + e~ — > 7*(p) — >• 
TT+(q 1 )TT-(q 2 )-i(k,e), see Fig. [TJ 





Fig. 1. Feynman diagrams of e + e 
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The amplitudes equivalent to those given in Ref. 12 1 
are obtained. It can be written as 



(1) 



where F M = v(pi, Ai)7 M u(p 2 , A 2 ). The p u X 1 ,p 2 , A 2 are 
momenta and helicities of the incoming electron and 
positron. V M defines the spin state of the intermediate 
!*■ 

Let us focus on the part for virtual photon decay. 
Following conventions of Ref. 12j, the final interac- 
tion part of the Born matrix element for such process 
is 

eF 2 „{p 2 ) 



HS{p,qi,q 2 ) 



P 



(2) 



Here p = q± + q 2 . If photon is present, this part of the 
amplitude reads: 



W = ^^l(q 1 + k-q 2 r^ 
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It makes sense to rewrite Eq. Q explicitly as sum of 
two gauge invariant terms: 

fP = JTj 4 + JT& 1 (4) 
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(6) 



P° \ q 2 -k + q 1 -k 

One can easily see that Eq.© has a typical form 
for amplitudes of QED exclusive exponentiation (l3 |. 
that is Born-like -expression multiplied by an eikonal 
factor {^f^ — ^t) - ^ e ex P ress i° n m front of the 
factor indeed approaches the Born one in soft pho- 
ton and collinear photon limit. Thus, it is consistent 
with LL level factorization into Born amplitude and 
eikonal factor. 

If one takes separation ((4]) for the calculation of 
two parts of spin amplitudes, then after spin average, 
the expression for the cross section takes the form: 
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We should stress that Eq.Q can have its first term 
even closer to Born-times-eikonal-factor form. For 
that purpose it is enough to adjust normalization 
of the first part of Eq.© to Born amplitude times 
eikonal factor, and replace |Mj| 2 with 



\M' I \ 2 = \M I Y 
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(8) 



Compensating adjustment to the remaining parts of 
Eq.© is then necessary. Since J2\ £ \ M'i\ 2 is the ex- 
pression used in PHOTOS Monte Carlo in Ref. such 
a modification is of interest. In the next section, we 
will perform our numerical investigations with respect 
to Ref. Q which is a reference for us. 

2.2 Numerical results 

We will show results at 2 GeV center of mass en- 
ergy. Comparison of result from ^ A f \M'j\ 2 with re- 
sult from PHOTOS with matrix element taken from Ref. 
!7i] is shown in Fig 121 One can see that agreement is 
excellent all over the phase space. It is true only for 
the case when distributions are averaged over the ori- 
entation of the whole event with respect to incoming 
beams (or spin state of the virtual photon) . 

If instead of Y^\ t \^-'i\ 2 one would use directly 
^ Ae |Mj| 2 , that is when normalization of Born-like 
factor is not performed, difference with respect to for- 
mulas in Ref. Q is much larger, see Fig. |3] 
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Fig. 2. Comparison of result using 
| | 2 (green line) with that using 
matrix element taken from Ref. [3| (red line). 
Black line represents their ratio. 




Fig. 3. Comparison of result using 
\Mi | 2 (green line) with that using 
matrix element taken from Ref. [3] (red line). 
Black line represents their ratio. 



be understood as bremsstrahlung, but rather as gen- 
uine process. Anyway in that region of phase space 
scalar QED is not expected to work well. Note that 
the difference between results of Figs. [5] and 0] is only 
0.2% of the total process rate. That is why our de- 
tailed discussion is not important for numerical con- 
clusions, but important for understanding the under- 
lying structure. 

From Fig. [5] one could conclude that the univer- 
sal kernel in Ref. [7], for arbitrary large samples, 
is equivalent to the matrix element J2\ c l-^rl 2 - But 
differences appear in distributions sensitive to initial 
state spin orientation, see Figs. [5] and H3 On these 
plots angular distributions of the photon momentum 
and 7r + momentum with respect to the beam line are 
shown, respectively. Regions of phase space giving 
near zero contribution at the Born level are becom- 
ing more populated if approximation for the photon 
radiation matrix element [7| is used. 




Finally let us compare result of complete scalar 
QED matrix element with that of matrix element 
taken from Ref. [7J, see Fig. At high photon energy 




Fig. 4. Comparison of result using complete 
matrix element (green line) with that using 
matrix element taken from Ref. Fjj (red line). 
Black line represents their ratio. 

region, there is clear surplus of events with respect 
to formula in Ref. Ml. That contribution should not 



Fig. 5. Comparison of result using 
£\ I A/j| 2 (green line) with that using 
matrix element taken from Ref. [?| (red line). 
Black line represents their ratio. 
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Fig. 6. Comparison of result using 
£\ e | A/} | 2 (green line) with that using 
matrix element taken from Ref. [3] (red line). 
Black line represents their ratio. 
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± — > 7T+7T e ± 



1/(7) 



In this section radiative corrections to A e4 decay 
are calculated following approximations explained in 
Ref. 14] . We compare it with Coulomb corrections 
used by NA48/2 collaboration. The hard photon 
bremsstruhlung is calculated analytically with soft 
photon and collinear photon approximation. Numer- 
ical tests with PHOTOS Monte Carlo are performed. 

3.1 QED radiative corrections to A e4 decay 

Consider A e4 decay, 

A± (p) -> 7T+ (q + ) + 7T- (q_ ) + e± (p e ) + i/(p„ ) . (9) 

In the framework of Scalar QED but neglecting dia- 
grams with photons emission from hadronic or weak 
blocks, one can calculate the virtual photon correc- 
tions. Note that the electron photon vertex is taken 
from standard QED. Contribution of virtual diagrams 
reads 





a.) b.) 
Fig. 7. example vertex diagrams [l^ 
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where to is the charged pion mass, A is photon mass 
used as infrared regulator, A is ultraviolet (UV) cut 
off. In EqJIOp we have defined 
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A„ depends on masses of particles and kinematics. 
Note that the UV-divergent part in Eq. (fT0|) will can- 
cel if the renormalized coupling G 2 F V 2 S is used instead 
of bare ones, 
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The soft photon contribution can be easily ob- 
tained, resulting in expression which is gauge invari- 
ant: 



dT s . 
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here k is photon momentum, u> is photon energy. It 
is straightforward to integrate out solid angle of pho- 
ton momentum k and over its energy w up to a limit 
lu < Ae; we obtain: 
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(14) 



Function A" s is dependent on masses of particles and 
kinematics. Soft singularity is again regularized with 
the photon mass A. 

The contribution of soft and virtual photons can 
be easily combined. It reads 



<£r B orn+virt+soft - . n . 7ra( 1 + /3 2 ) 

= 1 + aFs + ■ 



dT Bol . n 



+-K VS 

IT 



(15) 



where 



Ps = 21n| + | <r=^(2 P -l), (16) 

the expression of K V3 depends not only on masses 
of particles and kinematics, but also on soft photon 
energy cutoff Ae. 

Staring from a certain energy threshold photons 
can be observed, at least in principle. Such contri- 
bution is called hard photon radiation and we will 
assume that it is defined by condition uj > Ae. If 
soft photon matrix element Eq. (114[) is used, and max- 
imum energy of photon is taken from kinematical 
constraint E e , then Eq. (fl"4"|) can be used to calcu- 
late leading double logarithmic result of hard pho- 
ton radiation. For that reason one need to subtract 
Eq. (ll4]) calculated for E e from that for Ae. Thus the 
contribution of hard photon emission by soft-photon- 
approximation is obtained, 

= -In — U+ — - — t-2p 



Hard 3oft _ likc 
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2Pe-q+ 
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(17) 



This result has to be corrected for single logarithm re- 
lated to collinear photon emission along the charged 
outgoing decay products. 
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Photon emission from e will give collinear singu- 
larity. Remaining charged products, ir + and n~ are 
not relativistic. The electron part can be calculated 
with the help of collinear-photon-approximation. The 
squared collinear amplitude can be found in Ref. flj| | . 
The actual formula for collinear contribution is taken 
from Ref. \li 



r/r 



Colli] 



here z 



4tt 2 



d 3 k 



1 + z 2 



uj p e -k \1 — z p e -k 



> (18) 



— -f—. First we integrate over solid angle 
of photon direction, later we integrate photon energy 
from Ae to the maximum value E e . Then the result 
of hard photon emitting from e is obtained using 
collinear-photon-approximation, 



= -aP» 



(19) 



Eq. (fT9"|) includes term proportional to 2p— 1 which is 
already present in Eq. (jT7)) . This would lead to dou- 
ble counting. That is why we remove this term from 
Eq. (fT7|> and obtain 



IT \eJ 
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(20) 



Finally hard (real) photon bremsstrahlung for pho- 
tons of energy above Ae reads 



cT Ha 



ard c 



dr Ua 



(21) 



If one adds real and virtual photons contribution 
together, one can obtain the expression as following, 



lorn +virt + real 



1 



to(1 + /3 2 ) 
2(3 



7T 



(22) 



Complicated, but numerically small function K is de- 
pendent on masses of particles and kinematics of this 
process. Note that our final analytical result does not 
depend on the large logarithm In , or soft photon 
energy cut Ae. 

As one can see, Eqs. (|15[) (f2Tj) . (|22|) are obtained 
with the help of approximations. Effectively it was 
assumed that matrix element at Born level can be al- 
ways factorized out and photonic corrections can be 
calculated independently. Further corrections are as- 
sumed to be negligible and not affecting the substan- 
tial nature properties of hard interaction. This may 
be good as starting point, but cannot be left without 
future discussion/improvement£|. 



Our formulas are based on the same scheme of 
calculation as explained in Ref. [14| and in principle 
they should coincide numerically. Some differences in 
analytical results are still present. Also some numer- 
ical results still remain different. The exact expres- 
sions for K s , K v , K vs and K, as well as differences 
between our analytical results and these in Ref. [14 1 
will not be listed here for the limit of paper length. 
They will be present elsewhere. 

3.2 Numerical tests 



Let's start with Eq. (|2"2"|) . the correction for distri- 
bution when photon is integrated out but all other 
kinematic variables are kept. In Figs. |8] and [9] 
we show that dominant part of Eq. (|2"2"|) represents 
Coulomb correction. The difference is much smaller 




Fig. 8. Coulomb correction from Ref. [Io| . 
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Fig. 9. Radiative correction in Eq. (|22[). 



*We are grateful to Prof. J. Gasscr for stressing this Doint. 
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than the effect of Coufomb correction itself, see FigflOl 
where results for 1000 000 Born level events are placed 
in the histograms. We may conclude that our numer- 
ical implementation of Eq. ([^|) works well since its 
dominant part represents Coulomb correction. 
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Fig. 10. Difference between radiative correc- 
tion in Eq. (122 fl and Coulomb correction from 
Ref.[i"ol] calculated event by event. 

We have done numerical tests with PHOTOS for soft 
photons and found the distribution of soft photons 
from PHOTOS is as same as given in the soft photon 
expression fEq.([T3|'). We also have done the simi- 
lar test using hard photon expression (Eq. (|2Tj) ) and 
found it matches result from PHOTOS simulation in the 
soft photon region, too, as it is expected. 

For harder photons, we compare result of PHOTOS, 
with results of Eq. ([l4|) for soft photon and Eq. (|2Tj) 
for hard photon. Their ratios are given in Fig. 1111 
where effect as function of upper limit on photon en- 
ergy was used. 100 000 000 PHOTOS event samples 
were generated using single fixed Born level event, 
i.i . , , . 




0.5 1 ■ ■ 1 

0.0001 0.001 0.01 0.1 

E 2 /sqHS 

Fig. 11. Ratios of result from PHOTOS to results 
using soft photon Eq. (|14p and hard photon 
Eq. PT) . Here E\ < E 1 < E2 where we set 
-t=t- = 2.5, E2 varies from 0.0002V 7 S to 0.07^. 
Note that the energy of electron E e w 0.07\/S. 



We can conclude that agreement is good, as ex- 
pected. Though differences especially in harder pho- 
ton energy ranges can be seen. However even at the 
end of the spectrum, where distribution is poorly pop- 
ulated, differences are at 10 % level only. 

4 Summary 

We have presented the new tests of PHOTOS Monte 
Carlo, where the exact matrix element of 7* — > 7r + 7r~7 
is implemented and its numerical result is compared 
with the kernel of PHOTOS. QED radiative correction 
to process K ± — > l ± i/Tr + Tr~( r y) is also studied ana- 
lytically. Reasonable numerical agreement with sim- 
ulations including Coulomb correction and PHOTOS 
Monte Carlo was found. Since several assumptions 
are employed, further work is necessary. Our result 
is of practical interest for experiments. They con- 
firm that at least on technical level the Monte Carlo 
program works well as expected. 
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